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Abstract 


A simple 3 DOF analytical aerodynamic model of the Langley 
Winged-Cone Aerospace Plane concept is presented in a form suitable 
for simulation, trajectory optimization, guidance and control studies. 
This analytical model is especially suitable for methods based on vari- 
ational calculus. Analytical expressions are presented for lift, drag and 
pitching moment coefficients from subsonic to hypersonic Mach num- 
bers and angles of attack up to ±20 deg. This analytical model has 
break points at Mach numbers of 1.0, 1.4, 4.0 and 6.0. Across these 
Mach number break points, the lift, drag and pitching moment coeffi- 
cients are made continuous but their derivatives are not. There are no 
break points in angle of attack. The effect of control surface deflection 
is not considered. The present analytical model compares well with 
the APAS calculations and wind tunnel test data for most angles of 
attack and Mach numbers. 
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1 Nomenclature 


A 

C D 

Cr,e 

Cm 

C D ,c 

Cdj 

C Di 

C L 

Cl* 

Cm 

Cm 

Cn* 

CN,aa 

D 

ki,k2 
Kw(B) 
K b {W ) 
Kf 
L 

Lref 

Ln 

M 

T 

R 

Sb 

Sw 

Seff 

V B 

X 

/ 

X 

ac 

a 

P 

7 

A/ e 

A c /2 


aspect ratio of the wing 

drag coefficients 

root chord of the exposed wing 

base drag coefficient 

crossflow drag coefficient of the body 

skin friction drag coefficient 

induced drag coefficient 

lift coefficients 1/2 ^4^ 

__ dCj, 

pitching moment coefficient 

normal force coefficients N^j0F^rge 

_ dC N W 

da 

second derivative of normal force coefficient with respect to a 

aerodynamic drag 

apparent mass coefficients 

interference factor of wing on body 

interference factor of body on wing 

sonic leading edge correction factor 

aerodynamic lift 

reference length, equal to mean aerodynamic aerodynamic chord 

length of the conical nose 

Mach number 

local radius of the body 

maximum radius of the body 

maximum cross sectional area of the body 

reference wing area (planform area) 

flat plate area of the wing 

volume of the body 

axial distance measured from the leading edge of the body 
distance of aerodynamic center from the apex of exposed wing 
angle of attack 

CM 2 -T 

ratio of specific heats 

leading edge sweep angle of the wing 

mid chord sweep angle of the wing 
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p 


density of air 


Acronym 

APAS Aerodynamic Preliminary Analysis System 
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2 Introduction 


The National Aero-Space Plane (NASP) program envisages the develop- 
ment of a manned, single- stage-to-orbit vehicle capable of horizontal take- 
off and landing using conventional runways. Towards this objective, NASA 
Langley has proposed a simple winged-cone configuration (Fig.l) as a base- 
line model for simulation, trajectory optimization, guidance and control 
research. Refernce [1] presents a comprehensive tabular aerodynamic model 
using APAS [2]. However, for optimal control studies using variational cal- 
culus the tabular data is not very convenient. It is desirable to have an 
analytical aerodynamic model. 

The purpose of this report is to present a simple, analytical model of 
Langley winged-cone aerospace plane concept suitable for 3 DOF simulation 
and trajectory optimization studies especially based on variational meth- 
ods. Analytical expressions are presented for the estimation of lift, drag and 
pitching moment coefficients for subsonic, transonic, supersonic and hyper- 
sonic Mach numbers and angles of attack from 0 to ±20 deg. The analytical 
model presented here is based on DATCOM [3] methods and is supplimented 
by the theoritical methods wherever they are available. However, a formal 
difficulty arises because different estimation/ analytical methods have to be 
used at subsonic, supersonic and hypersonic speeds and this process intro- 
duces several break point Mach numbers. Across these break point Mach 
numbers, the aerodynamic coefficients become discontinuous. Using engi- 
neering approximations, the estimated aerodynamic coefficients have been 
made continuous across these break point Mach numbers. However, the 
derivatives of the aerodynamic coefficients with Mach number remain dis- 
continuous across these break point Mach numbers. The analytical model 
presented in this report has four break points at M — 1.0, 1.4, 4.0 and 6.0. 
There are no break points in the angle of attack. The effects of sideslip, roll 
angle and control surface deflection are not considered in this study. The 
Langley winged-cone configuration has a forward mounted canard which is 
supposed to be operative only at subsonic speeds and retracted at speeds 
exceeding Mach one. 

The present anlytical model is shown to be in fair agreement with the 
available wind tunnel test data [6] and APAS calculations [1]. This anlytical 
model may be used for 3 DOF simulations along with mass, inertia and 
propulsion models presented in [1]. 
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3 Estimation of Lift Coefficient 


In the following, we discuss the determination of normal force coefficient. 
Note that Cl = Cn cos a. 

3.1 Subsonic and Supersonic Speeds 

The normal force coefficient of the vehicle is given by, 

Cn = Cn,wb + Cn,c ( 1 ) 

where, Cn,wb and Cn,c are the normal force coefficients of the wing-body 
combination and canard respectively. The DATCOM [3] gives the following 
expression for the normal force coefficient of the wing body combination at 
subsonic and supersonic speeds: 

Cn,wb = Cn,n + [■ Kw{B ) + Kb{W)\ CN,e (2) 

where, Cn,n is the normal force coefficient of the nose (body), Kw(B ) and 
Kb(W ) are the interference factors of wing on body and body on wing 
respectively, Cn,€ is the normal force coefficient of the exposed wing, S e is 
the exposed wing area and Sw is the wing plan form area. For the given 
wing-body configuration, the exposed wing is defined as the parts of the 
wing outboard of the largest body diameter at the wing body intersection 
[3] (Figure 2) . From Table I, we have Sw = 3600 ft 2 and the exposed wing 
area S e was estimated as equal to 1176.8310 ft 2 . 

In the following, we discuss the determination of Cn.n^Cn^ li 'w(B), 
Kb(W ) and Cn,C at subsonic and supersonic speeds and match these pa- 
rameters at the break point Mach numbers. 

3.1.1 Estimation of Cn,n 

At subsonic speeds (M < 1.0), from DATCOM [3] 

C n ,n = 2 (fc 2 - ki)«|r (3) 

where the apparent mass coefficient, k 2 — ki, is equal to 0.92 for the 
subject configuration, Sb is the maximum cross sectional area of the body 
and a is the angle of attack in radian. In this report, unless otherwise stated, 
a is in radian. For the given configuration, we have Sb = 520.573 ft 2 . The 
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equation ( 3) is based on linear slender body theory and hence, we do not 
have nonlinear terms in the expression of Cn,N at subsonic speeds. 

At low supersonic speeds (1.40 < M < 4.0), the DATCOM [3] gives: 

Cn,n = 2(k 2 - ki)a-~- 4 — (4) 

ow jw 

where, the first term on the right hand side is the same linear slender 
body theory term of equation ( 3), Cd, c is the cross flow drag coefficient 
and S p is the planform area of the body which was estimated as equal to 
1893.1770ft 2 for the given configuration. It may be noted that the slender 
body result of equation ( 3) is supposed to be independent of Mach number. 
The nonlinear term on the right hand side of the equation ( 4) is due to the 
cross flow drag of the body. In general, the cross flow drag coefficient Cd, c 
depends on Reynolds number, Mach number and angle of attack. However, 
for simplicity we ignore the dependence of Cd, c on angle of attack and instead 
use a mean angle of attack of 10 deg to evaluate Cd, c at some selected Mach 
numbers in the range of 1.4 to 4.0. This data was then represented in the 
following analytical form: 

Cd,c = ciAf 3 -{■ c 2 M 2 4- c 3 M 4- C 4 (5) 

where C x = 0.16018212e - 01, C 2 = -0.21232216, C 3 = 0.83332039 and 
C 4 = 0.48085365. 

At transonic (1.0 < M < 1.40) and high supersonic Mach numbers 
(4.0 < M < 6.0), we use interpolation to obtain Cn,n • The interpolated 
transonic normal force coefficient is matched with subsonic value at M = 1.0 
and low supersonic value at M = 1.4. Similarly, the interpolated high 
supersonic normal force coefficient is matched with the low supersonic data 
at M — 4.0 and hypersonic data at M = 6.0. The interpolated expressions 
are: 

for 1.0 < M < 1.40, 

C NN = 1 - 8 ^ a5ij + 1.6766511 (M - \)o? (6) 

bw 

for 4.0 <M< 6.0, 

Cn,n = aia 3 4- a 2 a 2 4- a 3 a + a 4 (7) 

where, 

a x = 0.24197221M - 0.9678888 (8) 
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^2 

«3 

0*4 


— 0. 32024002 Af + 2.0393596e - 01 
— 0.55238185e - 01 M + 0.48702338 
1.5351068c - 04 M - 6.1404271e - 04 


( 9 ) 

(10) 

( 11 ) 


3.1.2 Estimation of CV,e 

For the exposed wing, at subsonic and snpersonic speeds, DATCOM [3] gives 
the following relation: 

C*N,e = h CjV.aa sin a | sin a I (12) 

where, C/va,e is the slope of normal force coefficient of the exposed wing 
with respect to a and Cn,uo is the second derivative of the normal force 
coefficient with respect to a. We assume Cjv a , e = Cz,a,e- 
For subsonic (0 < M < 1.0) Mach numbers, 



where, A is the wing aspect ratio equal to 4 tan A for a delta wing. Here, 
A — 90 - Ai e and (3 = y/M 2 - 1. Approximating the subject wing as a 
flat plate delta wing, we obtain k = 1. The parameter A c /2 is the sweep 
angle of the mid chord and is equal to 2 t * nA . Assuming A to be small so 
that tanA = A, which is justified for a highly swept delta wing like the one 
considered here, the equation ( 13) can be written as: 


47tA 


C Lo,,e x + 2 ^ A 2 (j _ M2) + 0 ,5 


(14) 


From Table I, we have Aj e = 76.0 deg. so that approximately, A = 0.25 
radian. Using 7r = 3.1428, we obtain: 

For 0 < M < 1.0, 


3.1428 

L “’ e_ 1 + 0 . 579 ^ 1 ^ 


(15) 


At transonic speeds (1.0 < M < 1.40), we use interpolation so that the 
interpolated coefficient matches with subsonic value of equation ( 15) at 
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M — 1.0 and supersonic value at M = 1.4 given by equation ( 17). This 
interpolation formula valid for 1.0 < M < 1.40 is given in the following: 

C Na , e = 0.35699257 M + 0.94479781 (16) 

At low supersonic speeds (1.4 < M < 4.0), using DATCOM [3] proce- 
dure, we obtain the linear term in normal force coefficient as follows: 


C Na , e = Ay(1.56250— 0.140620\/M 2 - 1) 
and for high supersonic speeds, 4.0 < M < 6.0, 


OjVa,« 


4.0 K f 
y/M 2 — 1 


(17) 


(18) 


In the above equations, Kf is the sonic leading edge correction factor. 
Using DATCOM [3] data and enforcing the matching of CjVa,e at M = 
4.0 and 6.0 and curve fitting the numerical data, we obtain the following 
expressions for Kf\ 
for 1.40 < M < 4.0, 


Kf = —0.0086 M 2 - 0.0385 M + 1.08470 (19) 


for 4.0 < M < 6.0, 

Kf = 0.25245975e - 01 M 2 - 0.15027659 M + 0.97881773 (20) 


Next we consider the evaluation of nonlinear term Cw^a- Using DATCOM 
[3] procedure, for subsonic speeds (0 < M < 1.0) we obtain, 

C N ,aa = 4.7783 - 1.4504CiVa fC (21) 


An approximation has to be introduced to evaluate Cj<r,aa at supersonic 
speeds because the given DATCOM [3] data do not cover the present config- 
uration for all combinations of Mach numbers and angle of attack of interest. 
The parameter Cjv,aa is almost constant with a at low supersonic Mach num- 
bers (1.0 < M < 2.0) but starts varying with a as Mach number increases 
above 2.0. At hypersonic speeds (as discussed later), this nonlinear term 
is proprtional to sina. Therefore, to retain this functional relationship, we 
will represent the available DATCOM [3] data as equations (22) and (23) 
with proper matching at M — 1.0, 4.0 and 6.0. (Note that the nonlinear 
term ( C^,aa ) doesn’t have a breakpoint at M = 1.4). 
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For 1.0 < M < 4.0, 

C Ni aa = -0.8467277 M + 3.736911 + 0.33333333 (M - 1) sin a (22) 
For 4.0 < M < 6.0, 

C N ,aa = -0.17500M + 1.05 + (1.90M - 6.6) sin a (23) 

3.1.3 Estimation of Kb(W) and Kw(B) 

From DATCOM [3], we find that for subsonic and transonic Mach numbers 
(0 < M < 1.4), Kw(B) = 1.35 and Kb(W) = 0.62. For supersonic speeds 
(1.4 < M < 6.0), using DATCOM [3] data and matching of the normal force 
coefficient of wing-body combination, we obtain the following expressions: 


K W (B) = -0.83998215c -01Af + 1.46759750c00 (24) 

jAb(PF) = q,\M^ + CL2M ^ -f- (13 + (I4M + 05 (2b) 

where, a x = 0.13778753e - 02, a 2 = -0.31571756e - 01, a 3 = 0.26725216e00, 
a 4 = — 9.9452836e — 01, and a 5 = 1.5698651e00. 

3.1.4 Estimation of Cn,c 

The canard is proposed to be deployed only at subsonic speeds to augment 
the static longitudinal stability and is supposed to be withdrawn immedi- 
ately for M > 1.0 [1], i.e., Cn,c = 0 for M > 1.0. Using DATCOM [3] data, 
for subsonic speeds (0 < M < 1.0), 

Cn,c = 0.08o-^- (26) 

bw 

where, angle of attack a is in deg and Sc is the canard area. From Table 
I, we have Sc = 154.3 ft 2 . Further, to mantain the continuity of Cn,c at 
M — 1.0, we gradually reduce the contributution of the canard to zero at 
M = 1.4 as follows, For 1.0 < M < 1.4, 


C N , c = (—0.20 M + 0.28 )a|2- 


(27) 
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3.2 Hypersonic Speeds 

Here, Cn,c — 0 so that Cn — Cn,wb ■ Further, at hypersonic speeds, the 
mutual interference between the wing and body can be ignored and the nor- 
mal force coefficient of the wing-body configuration can be assumed as equal 
to the sum of the contributions from isolated wing and body components, 
i.e., 

Cn,wb — Cn,b + Cn,w (28) 

According to DATCOM [3], the normal force coefficient of the body 
based on base area is given by, 

Cn ’ b = 7§ L KeH J d © (29) 

where, k is a constant and equal to 2 for pointed bodies, is the total 
body length, r is the local body radius, R is the maximum radius of the 
body, x is the local distance measured from the leading edge of the body 
and Ke(ot) is a function to be determined using DATCOM [3] data. 

To evaluate Ke(a), the integral in equation ( 29) was divided into three 
parts, conical nose, cylinderical center body and aft boattail. Using DAT- 
COM [3] data, the parameter I($ as a function of a was evaluated for each 
part and least square curve fitted to obtain analytical expressions. Then 
these contributions weere summed up to obtain the total Kd(&) for the 
body. The result of equation ( 29) was multiplied by to obtain the 
following expression for body normal force coefficient with respect to the 
reference wing area: 

Cn,b — «i« 3 + + a 4 (30) 


where, = 4.8394443e - 01, a 2 = 1.1791944e - 01, a 3 = 1.5559427e - 01, 
and a ,4 = 3.0702135e — 04. 

Reference [4] gives the following expression for the normal force of the 


wing based on flat plate approximation and applicable for small angles of 


attack: 


Cl,w — 2 a 2 


'2 (7 + l)Ma 

m Ma 6 


(31) 


Where 7 is the ratio of specific heats and for air, 7 = 1.4. We use equation 
( 31) to approximate the normal force coefficient of the wing as follows: 


‘4 sin a cos a . o 1 S e ff 

C N ,w= 77 + 0.8-M sin 3 a 

L M J Sw 


11 



Here, S e f / is the actual flat plate area of the wing and was estimated equal 
to 1357.6019 ft 2 . We have S w = 3600 ft 2 . Then, 

Cn,wb = Cn,b + Cn,w (33) 


4 Estimation of Pitching Moment Coeffficient 

The estimation of pitching moment coefficient at subsonic and supersonic 
speeds is based on the linear aerodynamic theory [3] and the concept of 

aerodynamic center . Once the aerodynamic center is determined, 

the pitching moment with respect to any desired moment reference point 
can be obtained. The moment reference center used in this report is shown 
in Fig 1. At hypersonic speeds, the method is more general and considers 
nonlinear variation of pitching moment coefficient with angle of attack. 


4.1 Subsonic and Supersonic Speeds 

In this section, we discuss the determination of C m at subsonic ( M < 1.0) 
and low supersonic (1.4 < M < 4.0) speeds using DATCOM [3]. The 
expressions for pitching moment coefficient at transonic (1.0 < M < 1.4) 
and low supersonic (4.0 < M < 6.0) speeds are obtained using interpolation. 
According to DATCOM [3], the location of the aerodynamic center of the 
wing body configuration is given by, 



where, 


ac 

c re 


N 


Cl«,n+ (Sf) , + ( 


J C L a,B{W) 
1 B{W) 


CLa,WB 


Clq,WB = C L a,N + Cl 0 ,W(B) + C La,B(W ) 


(34) 

(35) 


Here, J terms are the chordwise distances measured in exposed wing 

root chords from the apex of the exposed wing to the aerodynamic center 
and positive aft, Clcx,n = Cno,n » C Loii w(B) = Kw(B)CLa,e and C ia y W(B) = 

Kw(B)CLa,e • The determination of each of the is discussed in the 

following. 
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From subsonic to low supersonic speeds (0 < M < 4.0), the aerodynamic 
center location of the fore body as a function of the root chord of the exposed 
wing and referred to the exposed wing apex as the origin is given by, 



• -vr • 

CreC La,B 


(36) 


where, 

_ 2(^2 — &l) f^ N dS x ( f 0’~7\ 

Cma,B — y J X)dx ( 37 ) 

Here, Vb is the volume of the forebody, CLa,B = CNa,B — CNa,N, Ln is the 
length of the conical nose, x m is the distance of the moment reference point 
from the leading edge and S x is the local cross sectional area. As stated in 
DATCOM, we use x m — Ln • Evaluating the integral of the equation ( 37), 
we obtain, 
for 0 < M < 4.0, 

f^ac 

\ Cre 

For the given configuration, = 147.10 ft and the root chord of the 
exposed wing (c re ) was estimated as equal to 68.70 ft. The negative sign 
indicates that the aerodynamic center of the forebody lies ahead of the apex 
of the exposed wing root chord. 

Using DATCOM [3], the aerodynamic center of the exposed wing in 
the presence of body and as a fraction of the exposed wing root chord was 
estimated as, 
for 0 < M < 1.0, 

= —0.02915/? + 0.65 (39) 

W(B) 


where (3 = \/M 2 — 1. 

For 1.4 < M < 4.0, 

(£££] = 0.6650 (40) 

V Cre J W(B) 

The calculation of the aerodynamic center of the wing-lift carry over on the 

t 

body (f a£ )vr(B) involves several intermediate steps. Following DATCOM [3] 



N 


Ln 
3 c re 


(38) 
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we obtain, 
for 0 < M < 1.0, 


P^£ = 0.50 - 0.02387 6/3 (41) 

V Cre / B(W) 

and for 1.4 < M < 4.0, 

(^) = ai M 2 + a 2 M + a 3 (42) 

V Cre / B(W) 

where, a x = -1.4300077e-01, a 2 = 7.6328853e-01 and a 3 = -9.9911864e- 

02 . 

The pitching moment coefficient of the wing-body combination with ref- 
erence to the moment reference point (Fig. 1) is given by, 


c. 


mWB ~ 


CLa ,WBa(x' ac ~ Xref) 
Lref 


(43) 


where, x re j is the distance of the moment reference point from the lead- 
ing edge of the body and L re f is the reference length assumed equal to 
mean aerodynamic chord length. From Table I, we have, x re f = 124.0 ft 
and L re f = 80.0 ft. 

The contribution of the canard to the pitching moment coefficient at 
subsonic speeds (0 < M < 1.0) is given by, 


C m ,c — Cn,c 



(44) 


For M > 1.4, C mt c = 0. 

Then the total pitching moment coefficient of the vehicle at subsonic 
(0 < M < 1.0) and low supersonic (1.4 < M < 4.0) speeds is given by, 

C m — T C m ,c (45) 

At transonic (1.0 < M < 1.40) and high supersonic (4.0 < M < 6.0) speeds, 
we use interpolation. The transonic pitching moment coefficient matches 
with subsonic value at M = 1.0 and supersonic value at M = 1.4. Similarly, 
the high supersonic (4.0 < M < 6.0) pitching moment coefficient matches 
with low supersonic value at M = 4.0 and hypersonic value at M — 6.0. 
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For 1.0 < M < 1.40, 


C m = {a\M + a 2 )a (46) 

where, at = -7.12892780e - 01 and a 2 = 7.391297636c- 01. 

For 4.0 < M < 6.0, 


Cm = «1« 3 + + o 3 a + a 4 (47) 

where, 

al = — 3.8724006205d - 01M + 1.54896024820d0 (48) 

a2 = 2.2244718000d — 02M — 8.89788720000d — 02 (49) 

a3 = — 3.3111749650d - 02M + 1. 54278592200c* - 01 (50) 

o4 = - 1.5450500000c* - 06M + 6.18020000000d - 06 (51) 


4.2 Hypersonic Speeds 

As said earlier, at hypersonic speeds (M > 6.0), the mutual interference 
between the wing and body can be ignored and the pitching moment coeffi- 
cient of the wing-body configuration can be assumed as equal to the sum of 
the contributions from isolated wing and body components [4]. We divide 
the body into three parts, conical nose, cylinderical center body and aft boat 
tail. Then, 

s-i CN,N(Xref %cp,N^ . CiV,cy/(^re/ — cp,cyl ) 

- lZ 7 + 

G N,aft^ref %cp,aft) 

Lref 

where, Cn,n > CN tCy u and CN,aft are the normal force coefficients of the fore- 
body, cylinderical center body and the aft boat tail respectively. These co- 
efficients were evaluated using the equation ( 29). Further, x cp ^ = 2Xyv/3, 
%cp,cyi — T n T -^cyi/2 and 2Jcp,a/t — Lpj T Lpyi T L a ft/2. For the given 
configuration, Ln — 147.10 ft, Lcyl = 12.88ft and L a f t = 40.0ft. As be- 
fore, x re f — 124.0 ft. Evaluating equation (52), we obtain the following 
expression: 

C m ,B = ai« 3 + a 2 a 2 + a 3 a + a 4 (53) 
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where, ai = -3.9039436e - 02, a 2 == 4.4489436e - 02, a 3 = 5.7752634e - 02 
and a 4 — — 3.0901224e — 06. 

The pitching moment coefficient of the wing is given by, 


a 


m,W = 


l,w(,%ref * c cp,W’) 


' ref 


(54) 


where, Cn,w is given by equation (32) and x cPf w wa s estimated as equal to 
156.50285.’ Then, 

Cm ~ C m ,B T C m ,W (55) 


5 Estimation of Drag Coefficient 

The total drag coefficient of the vehicle is given by, 

Cd = Cnf y B 4 - Cof,w + Cpfy + Cob + Cdw,b + Cdw,w + 

Cdi,b + CDi,w + Cni,c (b6) 

where, Cdj,b is the skin friction drag coefficient of the body, Cdj,w is the 
skin friction drag coefficient of the wing, Cof,v is the skin friction drag co- 
efficient of the vertical tail, Cob is the base drag coefficient, Cdw,B is the 
wave drag coefficient of the body, Cdw,w is the wave drag coefficient of the 
wing, Cpi'B is the induced drag coefficient of the body, Coi,w is the in- 
duced drag coefficient of the wing and Cdi,c is the induced drag coefficient 
of the canard. We ignore the skin friction of the canard. For drag estima- 
tion purposes, we ignore the mutual interference effects. In the following, 
we discuss the evaluation of each of these components of the drag coef- 
ficients at subsonic, supersonic and hypersonic Mach numbers and match 
these components of drag coefficients at the break point Mach numbers. 


5.1 Body Skin Friction Drag Cofficient 

For subsonic speeds (0 < M < 1.0) DATCOM [3] gives, 


Cdj,b = Cj,b 


: 1+ (^ + 0 - 002 4 


A 

Sw 


( 57 ) 


where Cj^b is the skin friction coefficient of the body, l -f is the fineness 
ratio and S s is the wetted area of the body. For the given configuration, 
Ib = 200.0ft, d = 74ft and the wetted area ( S s ) was estimated as equal to 
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9426.7751 ft 2 . Cjb depends on body Reynolds number and Mach number. 
For simplicity, we use a mean Reynolds number (based on length) of 0.8xl0 9 
and further assume that the flow is fully turbulent over the entire body 
surface. Using DATCOM [3], we obtain C/b = 0.0015. With this, the 
equation (57) reduces to: 

Cd},b = 0.0044922 (58) 

In a similar fashion, Cof,B was evaluated at supersonic and hypersonic 
speeds. It was observed Cof,B at these speeds did not differ much from 
it’s subsonic value. Therefore, for simplicity, we use equation ( 58) for su- 
personic and hypersonic Mach numbers. 


5.2 Wing Skin Friction Drag Cofficient 

Using DATCOM [3], 


Cdj,w — Cf t w 


1 + L- + 100(-) 4 
c c 


Rls 


Jwet 

Sw 


(59) 


Where, Cf t w is the skin friction coefficient of the wing, ~ c is the thickness 
ratio, S we t is the wetted area and Rls is a parameter to be obtained using 
DATCOM [3]. From Table I, ~ = 0.04. We assume = 1. From DAT- 
COM [3], we get L = 1.2 for the given wing for which maximum thickness is 
located at 50 percent chord. For a mean wing Reynolds number of 0.4xl0 9 
based on root chord, we obtain Cf,w — 0.0022. Then, the equation ( 59) 
reduces to 

C D f,w = 0.0015 (60) 

As mentioned above, we also use equation ( 60) for supersonic and hypersonic 
Mach numbers. 


5.3 


Vertical Tail Skin Friction Drag Cofficient 


C Df,V = Cf t v 


Sv_ 

Sw 


(61) 


Assuming Cjy = 0.0022, we obtain 


C DftV = 0.00039 


(62) 


As in the case of the body and wing, we also use equation ( 62) for all 
supersonic and hypersonic Mach numbers. 
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5.4 Base Drag Coefficient 


From [5], for 0 < M < 1.0, 

C Db = [0.139 + 0.419(M - 0.161) 2 ] (63) 


For the given configuration, Sbase = 134.150 ft 2 . For supersonic and hyper- 
sonic speeds (M > 1.4), 



0.570] (Sbase\ 
M 4 \ V S w ) 


(64) 


For transonic Mach numbers (1.0 < M < 1.4), we use interpolation as 
follows, 

Cob = ciM + c 2 (65) 


where, cj = -1.80286235c - 01 and c 2 = 6.14229134e - 01. 


5.5 Body Wave Drag Coefficient 

For subsonic speeds, 0 < M <1.0, the body wave drag can be ignored. How- 
ever, to avoid the discontinuity in body wave drag coefficient at M = 1.0, 
we assume, 


Cdw,b — Cdw,b\M 6 ( 66 ) 

where, Cdw,bi is value of body wave drag coefficient at M = 1.0 given by 
the equation ( 67). At supersonic speeds, using DATCOM [3] the wave drag 
coefficient of the body was evaluated for 1.0 < M < 4.0 and the data was 
least square curve fitted to obtain the following expression: 


C D w,b = 0 2 (8.0101 


2.431M + 0.2443M 2 ) 


Sb_ 

Sw 


(67) 


From Table I, 9 = 5.0 deg or 6 = 0.08726 rad. As before, S B = 520.573 ft 2 
and Sw = 3600.0 ft 2 . For hypersonic speeds (M > 6.0), 

CW,B = 2 sinH (|5-) (68) 

For 4.0<Af <6.0, we use interpolation as follows, 

Cd»,b = (cxM + c 2 )|S- (69) 

bw 


where, c\ — -7.61321184e - 04 and c 2 = 1.97579419e - 02 and 
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5.6 Wing Wave Drag Coefficient 

For M < 0.60, Cdw,w — 0. Using DATCOM [3] and least square curve 
fitting the data, we obtain the wave drag coefficient of the wing for 0.6 < 
M < /- 05 a ,or approximately, 0.6 < M < 2.0, 

V COS 


- 0.000585A(M - 0.6) 

CDw ' W = 1.05 - 0.6VA 


(70) 


For the given configuration with A = 0.250, the equation ( 70) reduces to, 


C D w,w = 0.0018745(M - 0.6) (71) 

For y 1-05 . < M < — V- or approximately, 2.0 < M < 4.0, we have, 

VCOSAje — — COSAje 


C Dw,w = k cot Ale 



(72) 


For the given wing (^ = 0.04), the equation ( 72) reduces to, 

Cdw,w = 0.0016 (73) 


For — K— < M < 6.0, or approximately, 4.0 < M < 6.0, 

COS A(g 


„ 0.0064 

Dw ' w ~ Vm ^ I 


(74) 


However, a formal difficulty arises in using equations ( 71) and ( 74) 
because the wing wave drag coefficient has new break points at M = 0.6 
and M = 2.0 and the values of Cdw,W at M — 4.0 given by the equations 
( 73) and ( 74) are not identical. To overcome this difficulty, we approximate 
Cdw, w as follows so that we have the break points at M = 1.0, 1.4, 4.0 and 
6.0 as before: 
for 0 < M < 1.0, 


Cdw,w = Cdw,wiM ( 


(75) 


where, Cdw,Wi is the value of wing wave drag coefficient at M = 1.0 given 
by the equation (71). 
for 1.0 < M < 1.4, 


Cdw,w = ci M + c 2 


(76) 
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where, Cl = 2.12550000e - 03, and c 2 = -1.375700c - 03. 
for 1.4 < M < 4.0, 

Cdw,w — 0.0016 

for 4.0 < M < 6.0, 


Cdw,w = C\M + C2 


where, c\ = 4.8e — 4, and c 2 = 3.52e — 03. 
for M > 6.0, 


Cdu>,w = 0.0064 


( 77 ) 


(78) 


(79) 


5.7 Aft Body Drag Coefficient 

For subsonic Mach numbers, 0 < M < 1.0, C^^ft = 0. Using DATCOM [3] 
and curve fitting the data, we have for (1.0 < M < 6.0), 

C D , a ft = -0.0002AT 3 + 0.003 M 2 - 0.016 M + 0.0342 (80) 

For M > 6.0, Co,aft is negligibly small and is ignored. However, to have a 
continuous variation of Co,aft break point Mach numbers of 1.0 and 6.0, 
we approximate the variation of Co,aft at subsonic and hypersonic speeds 
as follows: 


Cc.a/t = C d , o/ „M 4 ;0< M < 1.0 (81) 

= CD, o/t2 e‘ (M ‘ 6 ' 0) ;M> 6.0 (82) 

where, Cn,afn and Co,aft2 are the values of CD,aft given by the equation 
( 80) at M = 1 and M = 6.0 respectively. 

5.8 Body Induced Drag Coefficient 

The induced drag coefficient of the body at subsonic, supersonic and hyper- 
sonic Mach numbers is given by, 

Cdi,b — Cn,b<* (83) 


Note that C/v.js = Cn,n • 
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5.9 Wing Induced Drag Coefficient 

The estimation of the induced drag coefficient of the wing is based on linear 
aerodynamic theory. Note that Cl,w = CLa,eOt for subsonic and supersonic 
speeds and Cl,w = Cjy.w cos a for hypersonic speeds. From DATCOM [3], 
the induced drag coefficient of the wing for 0 < M < 0.8 is given by, 

Coi,w = L ’ W ( 84 ) 

ne av 

where, e av — 0.450. For 0.8 < M < 1.4, using DATCOM [3], 


Cm,w = 0.4C£ )W 

For 1.4 < M < 6.0 from DATCOM [3], 

0.4010 + VM 2 - 1 


Cm,w = 


4.1887 


C 2 

° L,W 


For M > 6.0, 


Cm,w = Cl, 


$eff 

Sw 


(85) 


( 86 ) 


(87) 


Once again a formal difficulty arises here because we have a new break point 
at M = 0.8. To avoid this, the data is approximated as follows, 


Coi,w = 


‘(1 - M 4 ) 


ire „ 


+ 0 AM 4 


Cl > w ] 0 < M < 1*0 


(0.35699257M+ 0.94479781)C£ >w q 1.0 < Af < 1.4 
C\ w \ 1.4 < M < 4.0 


0.4010 + \/M 2 - 1 


4.1887 


clCdiwa + bCDi W 6', 4.0 < M < 6.0 

6.0 

bw 


( 88 ) 

(89) 

(90) 

(91) 

(92) 


where a — —0.5 M + 3.0 and b = 0.5 M — 2.0. CoiwA and Cdiwg are the 
values of Cniw at M — 4 and M = 6 obtained using equation ( 90) and ( 92) 
respectively. 
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( 93 ) 


(94) 


For M > 1.4, 

Cm,c = 0 

6 Results and Discussion 

The pictorial, three dimensional plots showing the variation of lift, drag 
and pitching moment coefficients with Mach number and angle of attack are 
shown in Figures 3 to 5. From these plots, it is observed that the coeffcients 
vary smoothly with Mach number and angle of attack. Typical variations 
of base drag, wave drag and induced drag coefficients with Mach number at 
an angle of attack of 4.0 deg. are shown in Figures 6 to 8 to illustrate the 
continuity of these coefficients across the break point Mach numbers. The 
variation of total drag coefficient at a = 4.0 deg. is shown in Figure 9 along 
with APAS [1] and wind tunnel data [6]. It is observed that the present 
anlytical model compares well with both data. 

The variation of normal force coefficients of the body, wing and canard 
are shown in Figures 10 to 12 to illustrate the continuous variation of these 
coefficients across break point Mach numbers. The total lift coefficient at 
a = 4.0 compares well with APAS [1] and wind tunnel data [6] as shown 
in Figure 13. In comparison with APAS predictions [1], the lift coefficient 
predicted by the present analytical model is closer to the wind tunnel data 
[6] especially at high Mach numbers. The total pitching moment coefficient 
at a = 4.0 predicted by the present analytical model is shown in Figure 14 
and is compared with APAS [1] and wind tunnel test data [6]. It is observed 
that significant differences exist between these data. 

To illustrate the nature of comparson at higher angles of attack, the 
variation of total drag, lift and pitching moment coefficients at a = 12.0 


(95) 


5.10 Canard Induced Drag Coefficient 

For 0 < M < 1.0, 

Cdi,c = 0.909 Cn,c^^~ 
bw 

For 1.0 < M < 1.4, 

Cdi,c = (c\M + c 2 ) -Tr- 
ow 

where, c\ = — 1.35988656e — 02 and = 1.90384119e — 02. 
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deg. are shown in Figures 15 to 17 along with APAS [1] and wind tunnel 
data [6]. As before, the drag and lift coefficients compare fairly well and 
significant differences exist in pitching moment coefficient. In view of this, 
some degree of care and caution may be excercised in using the pitching 
moment coefficient data. 

7 Concluding Remarks 

We have presented a simple 3 DOF analytical aerodynamic model of the 
Langley Winged-Cone Aerospace Plane concept suitable for simulation, tra- 
jectory optimization, guidance and control studies especially for methods 
based on variational calculus. The present analytical aerodynamic model 
may be used along with the mass, inertia and propulsion models of [1]. 
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Table I. Geometric Characterstics of The Configuration 


Wing 

Reference (planform) area, ft 2 ... 3600.0 

Aspect ratio ... 1.0 

Span, ft ... 60.0 

Leading edge sweep angle, deg ... 76.0 

Trailing edge sweep angle, deg ... 0.0 

Mean aerodynamic chord, ft ... 80.0 

Airfoil section . . . diamond 

Airfoil thickness to chord ratio, (£) ... 0.04 

Incidence angle ... 0 

Dihedral ... 0 

Wing Flap 

Area each, ft 2 92.30 

Chord (constant), ft 7.22 

Inboard section span location, ft . . . 15.00 

Outboard section span location, ft . . . 27.78 

Vertical Tail 

Exposed Area each, ft 2 645.70 

Span, ft ... 32.48 

Leading edge sweep angle, deg ... 70.0 

Traibng edge sweep angle, deg ... 38.13 

Airfoil section . . . diamond 

Airfoil thickness to chord ratio, (|) ... 0.04 

Rudder 

Area, ft 2 ... 161.4 

Span, ft ... 22.80 

Chord to vertical tail chord ratio ... 0.25 
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Table I. Concluded 


Canard 


Exposed Area, ft 2 

154.30 

Aspect ratio 

5.48 

Span, ft 

33.6 

Leading edge sweep angle, deg 

16.0 

Trailing edge sweep angle, deg 

0.0 

Airfoil section 

NACA 65A006 

Incidence angle 

0 

Dihedral 

0 

Axisymmetric Fuselage 

Length 

200.0 

Cone half angle, deg 

5.0 

Cylinder radius (maximum), ft 

12.87 

Cylinder length ,ft 

12.88 

Boattail half angle, deg 

9.0 

Boattail length, ft 

40.0 

Moment reference center, ft 

124.0 
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Figure 1 Configuration geometry 
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Figure 2. Schematic illustration of the 
concept of exposed wing area 
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Figure 6. Variation of base drag coefficient with 
Mach number at a=4.0deg. 
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Figure 8. Variation of induced drag coefficient with 
Mach number at a=4.0 deg. 
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O APAS [1] 

Present Analytical Model 

O Wind Tunnel Data [6] 


\m 


Mach Number 


Figure 9. Variation of total drag coefficient with 
Mach number at a=4.0 deg. 











Mach Number 


Figure 13. Variation of total lift coefficient 
with Mach number at a=4.0 deg. 







Analytical Model 

O APAS [1] 

O Wind Tunnel Data [6] 
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Figure 14. Variation of total pitching moment coefficient 
with Mach number at a=4.0 deg. 
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Figure 15. Variation of total drag coefficient 
with Mach number at a= 12.0 deg. 
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Figure 16. Variation of total lift coefficient 

with Mach number at a=12.0 deg. 


42 




Mach Number 


Figure 17. Variation of total pitching moment coefficient 
•with Mach number at o<=12.0 deg. 
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